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Abstract 

We present an exhaustive analysis of transport measurements performed in 
twinned YBa2Cus07 single crystals which stablishes that the vortex solid- 
liquid transition is first order when the magnetic field H is applied at an an- 
gle away from the direction of the twin planes. We show that the resistive 
transitions are hysteretic and the V-I curves are non-linear, displaying a char- 
acteristic s-shape at the melting line H m (T), which scales as e(9)H m (T,6). 
These features are gradually lost when the critical point H*(9) is approached. 
Above H*{9) the V-I characteristics show a linear response in the experimen- 
tally accessible V-I window, and the transition becomes reversible. Finally 
we show that the first order phase transition takes place between a highly 
correlated vortex liquid in the field direction and a solid state of unknown 
symmetry. As a consequence, the available data support the scenario for 
a vortex-line melting rather than a vortex sublimation as recently suggested 
[T.Sasagawa et al. PRL 80, 4297 (1998)]. 
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I. INTRODUCTION 



During the last years, vortex physics in high-temperature superconductors has become a 
major topic of research. The main reason for this interest is that the interplay between ther- 
mal fluctuations, anisotropy and disorder determine the existence of several vortex phases 
in the magnetic phase diagram of these materials. These phases are separated by differ- 
ent kinds of thermodynamic transitions. In particular, it is now well established that in 
clean YE^CusOy crystals the vortex solid transforms into a liquid through a first order 
phase transition! - !. This thermodynamic transition shows up in the transport properties 
as a sharp drop or "kink" in the resistivity, p(T), at the melting temperature T m (H, 0)03, 
and has an hysteretic behaviour both in temperature and field! which corresponds to a 
superheating of the solid phase!. 

When correlated disorder is present in the sample, as for example in twinned YBa2Cu307 
crystals, the transition transforms into second order! and occurs between a solid called 
Bose- Glass!! and an entangled vortex liquid. The Bose-Glass phase, existing only within a 
small angular region around the direction of the correlated defects (8 = 0), is characterized 
by a universal behaviour of the non-linear and linear resistivity with well defined critical 
exponent At the Bose- Glass transition temperature, T BG (H, 9), the presence of correlated 
defects introduces changes in the thermodynamic properties of the mixed state that go 
beyond the mass anisotropy approximation. In contrast to T m (H, 6) which smoothly follows 
the angular dependence given by the anisotropy, Tbg(H,9) shows a sharp cuspi'iaround 
9 = 0. Within this small angular region, p(T) has a smooth temperature dependence 
near the transition, while at larger angles a similar kink to that observed in untwinned 
crystals developsBB!. Using the Lindeman criterion and the scaling rules for anisotropic 
superconductors! Kwok and co-workersEI fitted the angular dependence of the temperature 
at which the kink occurred. This fit was interpreted as an indication of a recovery of the 
melting transition in the vortex system when the magnetic field was tilted away from the 
planar defect directionB. 
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Using this interpretation, Langan et al. showed in a recent paperS that the kink in 
p(T) is suppressed for fields larger than a certain field H*(6). This behaviour is similar to 
that observed in untwinned crystals where a critical end point for the melting line has been 
identified^, and gives support to the occurrence of a first order solid- liquid transition when 
the field is rotated away from the twin planes. 

An interesting finding in oriented twinned crystals in inclined magnetic fields was done 
by Morre et al J, who showed that the transport properties of the liquid state at angles well 
beyond the so-called depinning angle, 9d, are quite different from those observed in clean 
samples. The depinning angle concept is commonly used in the literature to indicate the 
angle beyond which the twin boundaries becomes ineffective as a correlating potential. Morre 
et al. showed, however, that the twin boundary potentials continue to induce vortex velocity 
correlation even for 9 » Their results demonstrated that the vortex liquid remains 
correlated in the field direction above the resistivity kink temperature. This behaviour is in 
sharp contrast to that observed in untwinned samples where the vortex velocity correlation 
in the field direction is lost at the melting temperatureS The important difference between 
the characteristics of the two liquids was used to cast doubts about the interpretation of the 
resistivity kink in twinned crystals as a manifestation of a first order transition!!. 

In order to investigate further the nature of the vortex solid-liquid transition in the 
presence of correlated defects at finite angles, we studied the vortex dynamics in YE^CusOy 
crystals with oriented twin boundaries. We performed exhaustive transport measurements 
as a function of magnetic field, H, and angle, 9, with respect to the twin boundary direction. 
Characterizing the resistivity kink near the transition by its width, we show that it remains 
almost constant as H is increased until a field H*{9) is reached. Below this field the V- 
I characteristics are non-linear and, what is more important, the resistive transitions are 
hysteretic. Moreover, our results indicate that the hysteresis corresponds to a superheating 
of the vortex lattice, in agreement with the results obtained by Charalambous et al.i Above 
H* which scales as e(9)H*(9), the kink width suddenly increases, the V-I characteristics 
become linear and the hysteresis disappears. We find that the resistivity kink height follows 



a universal behaviour when the reduced variable e(8)H is used, in contrast to the results of 
Langan et al.E_J. Our data together with those previously reporte d@a lead us to conclude 
that when H is tilted away from the direction of the twin planes the transition to the vortex 
liquid state, characterized by a steep jump in the resistivity, is indeed a first order phase 
transition. This transition occurs between a highly correlated vortex liquid! and a solid of 
unknown symmetry and therefore supports an interpretarion in therms of line-like melting)^! 
rather than a sublimation of the vortex latticeB. 



II. EXPERIMENTAL 

We carried out transport measurements on YBa 2 Cu307 single crystals with only one 
family of oriented twins. The crystals were grown using the self flux methodEl. The fully 
oxygenated samples had typical dimensions 1 x 0.3 x 0.01mm 3 , critical transition temper- 
atures T c = 93. 2K, and transition widths AT C < 0.6A". Four parallel contacts separated 
by 150/im were made with silver epoxy over evaporated gold pads, resulting in contact re- 
sistances lower than 1Q. The crystal was mounted onto a rotatable sample holder with 
an angular resolution better than 0.05° inside a commercial gas flow cryostat with an 18T 
magnet. The current was injected at 45° off the twin planes. Transport measurements 
were performed using conventional DC techniques. All resistivity measurements were made 
within the linear regime using a current density J < 3A/cm 2 . 



III. RESULTS AND DISCUSSION 

In Figure |] we show the angular dependence of the temperature at which the resistivity 
of the crystal becomes zero within experimental resolution^!, at an applied magnetic field of 
6T. Similar results were obtained for other applied fields. As clearly seen, there is a cusp at 
small angles which is indicative of the Bose-Glass phase, as has been recently demonstrated 
by Grigera et all. In this paper we will concentrate on the investigation of the nature of the 
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solid-liquid transition in the angular region away from that in which the Bose-Glass phase 
exists. 

In Figure ^| we present R(T) data for three different angles, as a function of the applied 
magnetic field. As can be seen, the transitions at low fields display a sharp drop to zero 
which we identify with the characteristic resistivity kink. Note that, for a fixed angle, the 
kink is washed out as the field is increased, an effect similar to that observed in ref. [10] 
and to that reported for untwinned crystals with the magnetic field applied parallel to the 
c-axisi. In the case under study, the field value at which the transitions starts to broaden 
depends on the angle between H and the correlated defects: the larger the angle, the larger 
the field at which the kink is washed out. This feature is more clearly seen in Figure |3| 
where, in order to quantify the above mentioned behaviour, we plotted the full width at half 
maximum (FWHM) of the temperature derivative, dp/dT, of the curves shown in Fig. |2| as 
a function of the applied field H. It can be seen that for 8 = 23° the transitions start to 
broaden at approximately H 12T, while this field value is increased up to 14T for 8 = 38°, 
and to almost 18T for 9 = 58°. The values of the magnetic field at which the transitions 
start to broaden are independent of the criterion used for the definition of the transition 
width. 

According to the scaling theory of Blatter et alJ, in an anisotropic superconductor a 
physical magnitude which is a function of angle and field should scale as the product He(9), 
where the anisotropy factor e{9) = ^ 'cos 2 (9) + / y~ 2 sin 2 (9), with 7 2 = m c /m a b. Such scaling 
for the FWHM, with 7 = 7, is shown in Figure [|, where we have also included data for other 
measured angles. Note that all curves collapse onto one and that the sudden increase of the 
transition widths occurs at the same reduced field H* = He (8) ~ 11T for all angles. The 
sudden change in the transition width is indicative of a corresponding change in the vortex 
dynamics below and above the characteristic field H*. 

As is widely accepted, one of the most powerful tools to investigate vortex dynamics is 
the measurement of V-I characteristics. We performed such measurements in a temperature 
interval around the transitions for magnetic fields below and above H*. Typical results for 
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the resistance as a function of the applied current at an angle 9 = 23° are shown in Figure 
|5]. In panel (a) we plot the measurements for a field of 8T, lower than the critical field 
H* (23°). At low currents and high temperatures the vortex response is ohmic, changing 
to a non-linear behaviour as the temperature is reduced towards T m = 80.4. Note that in 
the intermediate temperature region, the R(I) curves have a characteristic s-shape. This 
feature, which has been reported in clean samples where a first order phase transition in 
the vortex lattice was i dentifiedS, is characteristic of all the magnetic field region below 
H* (23°). It is more pronounced as H is reduced, becoming less evident as H approaches 
H*. Well above H* the R(T) curves are linear up to the maximum current that can be 
used without producing heating effects. In panel (b) we show the results for H ~ H* (23°). 
Clearly, the vortex response is remarkably different from that shown in panel (a). In a wide 
range of applied currents the response is linear. Non-linearities developed at high currents 
and low temperatures and are due to vortex loop excitations^. One may wonder if these 
non-linearities are a sign of a glass transition taking place at lower temperatures, a question 
that has also been raised concerning the nature of the transition above the critical point in 
clean crystals. If this were the case, the tails of the resistive transitions above the critical 
field should follow the scaling behaviour predicted by the theoryEl, R ~ (T — TgY^' 1 ^. The 
analysis of our resistivity data in terms of this scaling yielded negative results. We cannot 
discard, however, the occurrence of the glass transition, since this negative results might 
be related to the voltage resolution we have in our experiments. Safar et al.@ experimen- 
tally observed the above mentioned scaling in Bi 2 Sr2CaCu20 8+ 5 crystals by using SQUID 
picovoltimetry. 

The extensive experimental results shown above display special features pointing towards 
the occurrence of a first order solid to liquid transition in the vortex lattice when H < H* . 
In order to consolidate this scenario, we have performed measurements to look for hysteretic 
behaviour in the resistive transitions^. Due to experimental constrains, instead of searching 
for hysteresis in temperature, we performed such measurements at a fixed temperature 
(regulation better than 5 mK) and sweeping the magnetic field up and down. The results 

6 



for an angle 9 = 23° and a measuring current I = 50 /j, A are plotted in Figure || Panel (a) 
shows the data for T ~ T* where T* is the temperature corresponding to the critical field at 
9 = 23°. The arrows indicate the sense in which the field was sweeped. Within experimental 
resolution no hysteresis is seen in this region. However, when the temperature is increased 
in such a way that T > T* a clear hysteretic behaviour in the resistive transition develops, 
as can be seen in panel (b). It is important to mention that the width of the hysteresis is 
current independent between 50 and 100/iA When the applied current is increased above 
this value curve A shifts towards curve B and, at high enough currents (above 150/iA), the 
hysteresis is washed out. Following Charalambous et al.i we interpret this behaviour as 
indicative of a superheating of the vortex lattice. 

In Figure [F] we compare the superconductor phase diagram for our sample with that 
obtained in a clean untwinned crystal. The inset show the raw H-T data, for different 
angles and the data at zero angle for the untwinned crystal. In order to take into account 
the anisotropy change as the angle between H and the defects is changed, in the main panel 
we have used the corresponding scaling field He(9), while the reduced temperature scale 
t = T/T C (H = 0) is used to account for the different critical temperature at zero field of 
both samples. The collapse of all data on one universal curve not only provides an impressive 
graphical view of our interpretation of the resistivity kink temperature in twinned samples 
as a melting temperature, but indicates that the presence of the twin boundary potentials 
does not modify this temperature. Later on we will come back to this point, continuing now 
with the anisotropy dependence of other physical quantities. 

In clean crystals a relevant quantity related to the melting transition is the kink height 
at the melting temperature. It has been found to be angle and field independent^, and this 
particular behaviour has attracted theoretical interest^. In contrast, we have found that in 
45° oriented twinned crystals, the kink height is angle and field dependent, an observation 
already reported in ref. [H]. Since this quantity related to the occurrence of a first order 
phase transition which follows the anisotropy, we expect the kink height also to scale with 
it. In Fig. H we show such scaling. We have plotted the kink height measured at fixed 
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angles and increasing the magnetic field (open symbols), together with measurements of 
the same quantity but at a fixed applied field (6T) and reducing the angle 6*(full symbols). 
Clearly all data collapse onto a universal curve. The behaviour reported by Langan et al.0 
is distinctly different. The reasons for the lack of scaling in their data is related to the fact 
that their measurements at 9 = 5° were taken in an angular region were the dissipation at 
the transition is greatly reduced due to the effect of the twin planes on the vortex dynamics. 
This reduction is easily seen in V(9) measurements where a dip in the dissipation occurs&Bi! 
for 9 < 10°, an angle usually identified with 9d- One may wonder if this change in the vortex 
dynamics is a consequence of a corresponding change in the thermodynamic nature of the 
transition for 9 < 9^. This is certainly true when the Bose-Glass phase is reached since it 
has been demonstrated! that in this case the transition is second order. However, this phase 
existsB below 9bg ~ 2°, leaving a rather wide angular region (9bg < 9 < 9d) where the 
full kink height develops (see e.g. ref. H). One interesting possibility! is that this angular 
region comprises a reentrant liquid phase separating the Bose-Glass from the vortex solid, 
in a similar way to what has been predicted to occur near H cl in clean samples inil. In this 
case the resistive transitions below the developing kink should be linear due to an increased 
viscosity of the reentrant liquid phase as the temperature is lowered. Within this picture 9d 
would be the angle below which the reentrant phase starts to develope. Another possibility 
is that the angular region between 9bg an d 0d) is governed by critical fluctuations of the 
Bose-Glass phase associated to the new thermodynamic variable H± = Hsin{9). If this were 
the case the so-called depinning angle should be interpreted as the critical angle above which 
the first order solid-liquid transition sets in. Although a detailed analysis of the nature of 
the solid-liquid transition in this angular region is out of the scope of this paper, further 
investigation is underway to elucidate this interesting issue. 

In the following, we would like to comment on a few important points concerning the 
H — T phase diagram shown in Fig. [7]. The first one is related to the value of the critical 
field H* compared to that of the untwinned sample. It is well known that in untwinned 
crystals this field has not a universal value. Its magnitude depends on the amount of point 



like disorder present in each crystals, the larger the disorder the lower the critical field. 
The results in Figure ^ which show H* ~ Hg may therefore suggest that in tilted magnetic 
fields the presence of twin boundaries does not increase the amount of point like disorder. 
This suggestion seems to be corroborated by our measurements in another oriented twinned 
crystals and by the results shown in ref. with a value of H* ~ 11T, similar to that 
obtained in our crystals. 

The second remark is more fundamental since it is related to the nature of the solid and 
liquid vortex phases in the twinned crystal. In a recent paper Sasagawa et al.0 suggested 
that the first order solid-liquid phase transition in the vortex system of high temperature 
superconductors, including the less anisotropic YBa 2 Cu 3 07, corresponds to a vortex sub- 
limation rather than a line melting. Within this scenario the vortex velocity correlation 
length in the field direction should vanish at the melting temperature because the vortex 
liquid phase is formed by uncorrelated pancakes. As already mentioned in the introduc- 
tion, Morre et al. have showni, using the flux transformer contact configuration, that the 
vortex liquid in tilted magnetic fields mantains the vortex velocity correlation in the field 
direction even at rather large angles off the twin planes. Since we have concluded that the 
transition in tilted fields is indeed first order, their results imply that the solid transforms 
into a vortex liquid of correlated lines, contrary to what happens in clean samples. This 
has a two important implications: first the sublimation scenario proposed by Sasagawa et 
al.0 does not hold for the melting transition in twinned crystals for tilted magnetic fields, 
and second, although the twin boundary potentials do not affect the melting temperature 
they play an important role in building up the vortex velocity correlation in the field direc- 
tion. This last point indicates that the degree of correlation in the vortex liquid does not 
determine the nature of the thermodynamic transition. Therefore one may conclude that 
the liquid phase in clean and untwinned samples (in tilted magnetic fields) is essentially the 
same but with different dynamics due to the different vortex velocity correlation. On the 
other hand, it could also be possible that the nature of the liquid phase differs from that 
in untwinned samples. We speculate that for twinned crystals in tilted magnetic fields, the 
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entangled vortex liquid might be formed by stair-like correlated lines which are stabilized 
by the twin boundary potential. Within this picture, the solid might have a complicated 
ordered structure also formed by stair-like vortices. 

IV. CONCLUSION 

We have shown that the vortex solid-liquid transition in twinned YBa2Cus07 crystals 
for magnetic fields tilted away from the planar defect's direction is first order. This ther- 
modynamic transition is reflected in the transport properties as a sharp kink in the R(T) 
curves at the melting temperature T m (H,9). The resistive transitions are hysteretic below 
a critical field that scales with the anisotropy of the material, with the hysteresis indicating 
a superheating of the solid phase. The correlated vortex liquid phase near the melting tem- 
perature could be probably formed by stair like entangled lines that maintain their vortex 
velocity correlation in the field direction due to the effect of the twin boundary potential. 
The structure of the vortex solid is unknown. 
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FIGURES 

FIG. 1. Zero resistance temperature as a function of the angle 9 between the applied magnetic 
field (H = 6T) and the twin boundary direction 

FIG. 2. Resistivity as a function of temperature for different applied fields, (a) 9 = 23°, (b) 
9 = 38°, (c) 9 = 58°. The curves shown correspond to magnetics fields of H = 0, 2, 3, 5, 7,8, 9, 10, 
11, 12, 13, 14, 15, 16 and 18T. 

FIG. 3. Full Width at Half Maximum of the derivatives of the resistive transitions as a function 
of the magnetic field. 

FIG. 4. Full Width at Half Maximum of the derivatives of the resistive transitions as a function 
of the reduced magnetic field He{9). 

FIG. 5. V-I characteristics at an angle 9 = 23° for: (a) H = 8T, below the critical field, (b) 
H = 11T, just above the critical field. The V-I curves were taken at steps of 0.2 K. 

FIG. 6. Resistive transitions as a function of H at an angle 9 = 23°. Curves A were measured 
sweeping the field up, while curves B were recorded by sweeping the field down, (a) T = 78.52K, 
near the critical temperature, (b) T = 81.82K, above the critical temperature. 

FIG. 7. Comparison of the melting lines for the twinned crystal in tilted magnetic fields (open 
symbols) and untwinned crystal with H parallel to c axis after ref. [ll]. The main panel shows the 
scaled data, while in the inset in the inset the raw data are plotted. 

FIG. 8. Normalized height of the resistive kink at the melting temperature as a function of the 
reduced field. Open symbols correspond to measurement at fixed angles and varying the magnetic 
field. Solid symbols were measured at 6T varying the angle between 9 = 20° and 9 = 85°. 
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